function dev = std (x)
    N = size(x,'c');
    m = mean(x,'c');
    num = x - m;
    dev_ = sqrt(num * num' ./ (N-1));
    dev = real(dev_);
endfunction
